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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

THEORETICAL INVESTIGATION OF THE ATTACK PHASE OF AN 
AUTOMATIC INTERCEPTOR SYSTEM AT SUPERSONIC SPEEDS 
WITH PARTICULAR ATTENTION TO AERODYNAMIC AND 
DYNAMIC REPRESENTATION OF THE INTERCEPTOR 
By Windsor L. Sherman and Albert A. Schy 

SUMMARY 


This paper presents the results of a theoretical investigation of 
the attack phase of an automatically controlled interceptor. The 
specific objectives of this study were to investigate the dynamic repre- 
sentation of the airplane, the flight behavior of the airplane when it is 
a part of an interceptor system, and the effect of nonlinear aerodynamics 
on the airplane flight maneuvers. In order to conduct this study, a 
modern high-speed interceptor, a dynamically perfect radar-controlled 
director-type guidance system for a lead-collision course and a velocity- 
type automatic pilot were simulated on a large analog computer. Radar 
noise and gust effects were neglected. The results are presented as time 
histories of the airplane and system motions. 

The basic flight maneuver encountered is a 5 g climbing turn with 
roll and acceleration commands applied simultaneously. The use of a 
linear representation for the airplane caused errors in the airplane 
motion and, in addition, caused large errors to appear in the miss dis- 
tances. The principal effect of the nonlinear aerodynamics was to change 
the elevator and aileron motions. Secondary effects appeared in the 
longitudinal motions and sideslip and yawing velocities . 


INTRODUCTION 


The concept of automatically controlled interception of hostile 
aircraft was established when it became apparent that interception mis- 
sions would have to be carried out in all types of weather and that the 
operational problems of modern interceptor aircraft were making the per- 
formance of the interceptor mission increasingly difficult for the 
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■unassisted pilot. The philosophy of the automatic interceptor weapon 
system requires that the entire interceptor' mission, from take-off 
through landing, he accomplished under automatic control. Several mili- 
tary contracts for the development of operating interceptor systems have 
resulted in system stability studies in which the airplane is considered 
as a linear component of the system. None of the published studies, how- 
ever, have investigated the effects on system response of the nonlinear 
dynamic terms and nonlinear aerodynamic forces in the equations of motion 
of the airplane. The present study is primarily concerned with these 
effects . 

The attack phase was selected for this investigation because violent 
maneuvers may be required of the airplane. The effects of nonlinear terms 
in the airplane representation would naturally be most evident for such 
large motions. The objectives of the investigation were as follows: 

(a) To determine if significant errors are introduced in the air- 
plane motions and system response when a linear representation is assumed 
for an airplane and to determine the simplest set of nonlinear equations 
that can be used to represent the airplane in the interceptor system 
herein considered. 

(b) To study the effect of nonlinear aerodynamics on the airplane 
motions and system response. 

(c) To study the flight behavior of the airplane during an auto- 
matically controlled attack run for a system where gravity corrections 
are omitted from the roll command. 

In order to investigate the attack problem, a target, fire-control 
radar, fire-control and guidance computer, automatic pilot, and the air- 
plane were simulated. This required the use of a large analog computer, 
and through the cooperation of the U. S. Navy, the Typhoon Computer at 
the U. S. Naval Air Development Center, Johnsville, Pa., was made avail- 
able for this study. The use of this computing facility permitted the 
simulation of an advanced high-speed interceptor by the six-degree-of- 
freedom rigid-body equations of motion. In addition, aerodynamics that 
were nonlinear functions of the angle of attack and linear functions 
of Mach number were included in the airplane simulation. Radar noise 
and gust effects were neglected in this study. 

The attack run, which starts when the radar locks on the target and 
ends with the discharge of the interceptor's armament, was made by an 
interceptor flying at a Mach number of 2.2 against a target that had a 
Mach number of 1.4. This paper presents a description of the interceptor 
system assumed, results that illustrate typical interceptor motions, and 
the effects of changes in the dynamic and aerodynamic representation of 
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the interceptor on the motions of the airplane during the attack run. 

A summary of some of the more important results was previously presented 
in reference 1. 


SYMBOLS 


b 

c 


normal acceleration 
wing span 

mean aerodynamic chord 



E a> E e 

Wz 


g 

X X> X Y> I Z 


L 


M 

Ma,Me 


M* 

m 

N 

n 

Hrp 


unfiltered azimuth and elevation steering errors 
forces along interceptor body X-, Y-, and Z-axes 
acceleration of gravity 

moments of inertia about X, Y, and Z principal body axes 

K2 q constants 
rolling moment 

direction cosines between space axes and airplane principal 
body axes 

pitching moment 

azimuth and elevation miss distances in airplane principal 
body axes 

Mach number 

mass 

yawing moment 

interceptor normal acceleration, g units 
target normal acceleration, g units 
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a 

% 

R 

r 

S 

t 

t f 

t g 

u 

V f 

V m 

V s 


T 

v 

W 

w 

■X,Y,Z 

a 

(3 


airplane rolling velocity in principal body axes 
airplane pitching velocity in principal body axes 
local sonic dynamic pressure 

range 

range components in principal body axes 

airplane yawing velocity in principal body axes 

wing area 

time 

time of flight, interceptor 

time to go; time from present to firing point 

interceptor velocity along X body axis 
interceptor velocity 

missile velocity relative to interceptor 
velocity of sound, free stream 
target velocity 

interceptor velocity along Y body axis 
weight 

interceptor velocity along Z body axis 
right-hand cartesian coordinate system 
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angle of attack 
angle of sideslip 
flight-path angle 
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A 


when associated with variable, indicates perturbation 
in variable 


^a^e^r aileron, elevator, and rudder deflector angles, respectively 


e = V e a 2 + e e 2 

filtered azimuth and elevation steering errors in body axes 


€ a; e e 


®a> ®e 


fz 


SA 


*SE 


SR 


$ 

* 


g-limited flight-path command 

Euler angle (interceptor pitch) 
azimuth and elevation radar gimbal angles 

air density 

time of flight of rocket 
filter time constant 

aileron servomotor time constant 
elevator servomotor time constant 
rudder servomotor time constant 

Euler angle (interceptor bank) 

Euler angle (interceptor yaw) 

2 F X 


X-force coefficient, 


Y-force coefficient, 


pSVf* 


2 K 


pSV f 2 

Z-( normal) force coefficient, 


2Fr, 


pSV^ 
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rolling-moment coefficient. 


2L 

pSV f 2 b 


pitching-moment coefficient. 


2M 

pSV f 2 c 


yawing-moment coefficient. 


2N 

pSV f 2 b 


acx 

8& e 

dCy 

5T 



ac z 




8& r 




8C 


l 


2 V f 
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SC 


' n P 
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vPt> 


2V, 


Sc. 




n 


s^ 


2V 


Subscripts : 

b principal body-axes coordinate system 

c command, that is, $ c signifies roll command 

cr indicates critical value of variable 

i numerical subscript values 1, 2, 3 

l indicates limiting value of variable 

o indicates initial condition 

s ' indicates space coordinate system 

Dot indicates differentiation with respect to time. 

All angles measured in radians unless otherwise stated. 
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SIMULATION OF THE ATTACK PHASE OF THE AUTOMATIC 
INTERCEPTOR SYSTEM 


This study was made on the Typhoon Computer at the U. S. Naval Air 
Development Center , Johnsville, Pa. Reference 2 is a detailed descrip- 
tion of the Typhoon computer, whereas reference 3 discusses the simula- 
tion of this problem and presents a discussion of the interceptor system 
used. Appendix A of this paper summarizes the equations used to simulate 
the automatic interceptor attack problem for this study. 


Basic Assumptions 

The principal assumptions made in setting up this system are as 
follows : 

(a) The dynamics of the radar, computers, and flight data sensing 
instruments, such as rate gyros, are neglected. 

(b) Radar noise and gust effects are neglected. 

(c) Rocket ballistics are neglected. 

(d) Changes of atmospheric density with altitude are neglected, the 
value assumed being that for the altitude of the target as given by the 
NACA standard atmosphere (ref. k ) , 

(e) The angles a and p are small angles 

sin a ~ a sin |3 — p 

cos a ~ 1 cos p — 1 

(f) The angle e, the angle between the body reference system and 
principal axis, was small (e = 1 . 2 °) and was assumed zero. Thus, the 
body reference system was assumed to be coincident with the principal 
axes . 

(g) V f = \Jn 2 + v 2 + w 2 « — 

1 2 - a 2 

(h) Control surface servomotors are represented as first-order 
systems . 
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Components Simulated 

Figure 1 is a generalized block diagram of the interceptor system 
simulated for this investigation. The components of the system are tar- 
get, radar, fire-control computer, filters, command computer, autopilot, 
and airplane. The principal information flow routes have been indicated 
on the block diagram. 

The target .- The target was programmed to fly a straight-line course 
at an altitude of 50^000 feet. The target Mach number was l.k, which 
corresponds to an airspeed of 1,529 feet per second. Provisions were 
incorporated so that, if desired, the target could make a ±2g vertical 
plane maneuver that started at radar lock-on. 

The radar and fire-control computer .- The radar, which is mounted 
in the interceptor, tracks the target and supplies range, range rate, 
antenna gimbal angles, and the angular velocities of the line of sight 
to the fire-control computer. Since a perfect radar was assumed, the 
kinematic relations, equations (Al) to (At) of appendix A, were used to 
determine the required information. The geometry on which these kine- 
matics were based is presented in figure 2. The fire-control computer 
was a director-type computer with first-order prediction that solved the 
equations of a lead collision rocket-firing course for the miss distance 
ML NU , v 

parameters - — and - — Equations (A5) of appendix A were used 
§ ES 

to represent this computer. In deriving these equations it was assumed 
that the interceptor was armed with unguided rockets. For this study 
the rocket had an average velocity of 2,000 feet per second relative to 
the interceptor, and the time of rocket flight was 1.5 seconds. 

The command computer .- The miss-distance parameters, which were the 
output of the fire-control computer, were fed to the command computer. 

The command computer converts this miss-distance information to autopilot 
commands. The miss-distance parameters are first converted to steering 
errors by equations (A6), then filtered and corrected for cross-roll 
effects by equations (A 7 ) • These filtered steering errors are then con- 
verted to autopilot commands e e and 0 C by equations (A8) to (AlO) , 

where e Q is the normal-acceleration limited flight-path command and 0 C 

is the roll command. Both of these commands are uncompensated for the 
effects of gravity. 

The automatic pilot .- The automatic pilot was composed of two veloc- 
ity control loops, one for roll and one for the flight path. The roll con- 
trol system is described in reference 5 and the flight-path control sys- 
tem is described in reference 6. Figures 3 and k show the block diagrams 
of the lateral and flight-path control systems, respectively. The values 
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of the gains used are presented in appendix A as are the equations that 
describe figures 3 and k for simulation (eqs. (All) and (A12)). 

The roll control automatic pilot 'was the primary lateral control 
loop for the interceptor and was not changed during the course of the 
investigation. Two auxiliary regulators were considered for rudder con- 
trol. The first was a yaw rate feedback to the rudder for controlling 
the Dutch roll mode and the second was a p feedback that, in addition 
to providing Dutch roll damping, tends to reduce the sideslip during 
maneuvers . 

As indicated in reference 6, a low-gain longitudinal control system 
was satisfactory for tracking a nonmaneuvering target, whereas a high- 
gain control system was required to track a maneuvering target. In this 
study the low- and high-gain flight-path control system were used against 
nonmaneuvering targets in order to determine how the gain in the flight- 
path control system affected the airplane response when nonlinear aero- 
dynamics was included in the airplane simulation. 

In accordance with assumption (g), all servomotors in the automatic 
control system were represented by first-order equations, so that 

— £ = 1 with t q = 0.03 for all servomotors. 

&i 1 + t s D b 

The roll control system had an acceleration feedback of p. This 
feedback was included to alleviate control-surface rate-limiting oscilla- 
tions and is discussed in detail in reference 5 • 

Appendix B presents the modifications in the airplane effective 
damping and inertia characteristics introduced by the automatic pilot. 

The airplane .- The airplane was represented by the six-degree -of - 
freedom rigid-body equations of motion referenced to principal body axes 
(see eqs. (A13)). The terms in the brackets are the nonlinear cross- 
product terms normally neglected in the linear analysis of small motions. 
These terms were programmed so that they could be deleted from the equa- 
tions in order to study the effects introduced by representing the air- 
plane by linear equations. The interceptor was assumed to be in trimmed 
straight and level flight at a Mach number of 2.2, = 2,136 feet 

per second, at radar lock-on. This flight condition determined the 
initial conditions for the equations of motion and the direction cosine 
computation. The direction cosines are functions of the airplane angular 
velocities and were computed by equations (Al4) . The direction cosines 
were used to obtain the correct distribution of gravity forces in the air- 
plane equations and in the radar simulation for the coordinate transforma- 
tion between space axes and interceptor body axes. 
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The right-hand sides of the equations of motions (eqs. (A13)) are the 
forces and moments acting on the airplane. Equations (A15) and (Al 6 ) are 
the control- surface forces and moments, and equations (A17) and (Al 8 ) are 
the other aerodynamic forces and moments. The parameters a, angle of 
attack, 0 , angle of sideslip, and M*, Mach number, required in the com- 
putation of the forces and moments were computed by equations (A19) . All 
stability derivatives were functions of Mach number. The pitching -moment 
coefficient Cm and the stability derivatives CZp, Cnp, CZpj and ^ n p 

were nonlinear functions of the angle of attack in addition to being func- 
tions of Mach number. The stability derivatives were computed by equa- 
tions (A20) and (A21) , using servo multipliers. 


The analog for the airplane was prepared so that any of the following 
forms of the stability derivatives could be used: 

(a) constant derivatives 

(b) derivatives that vary with Mach number 

(c) derivatives that vary nonlinearly with angle of attack 

(d) derivatives that vary with Mach number and nonlinearly with 
angle of attack. 


The nonlinear angle -of -attack variations in C m , C^ , C 1 , C nR , 

p p P 

and C n p could be added to the problem singly or in any combination. 

As the estimates of the stability derivatives showed that was very 

small and, when compared with Cm^ and C^, would have practically no 


influence on the pitching -moment equation, it was omitted from the simu- 
lation. In the case of Cy^ Cz£, anc ^ no es tbnates or wind- 

tunnel results were available; therefore, these stability derivatives 


were neglected. Aerodynamic cross-coupling terms such as 
omitted from the simulation. 


HI 


P 


were 


Initial Conditions 

Five initial conditions were used in this study. (See fig. 5( a )-) 

Of these five initial conditions, two (i and II) were bow attacks and 
three (ill, IV, and V) were beam attacks. All of these initial conditions 
except V result in forward hemisphere attacks on the target. In initial 
condition V the interceptor starts as a beam attack that degenerates into 
a tail chase that came in 30° off the tail of the target. Figure 5( a ) 
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presents a sketch which shows the variables that define the initial con- 
ditions. The specific values of these parameters are given in the table 
of figure 5(a). Figure 5(b) is presented to clarify the lock-on orienta- 
tion of the target and interceptor for each initial condition. This figure 
shows the lock-on condition as it would be seen by an observer in an air- 
plane that is above and on the right rear quarter of the target. 

In these attack runs it was assumed that the radar had been tracking 
the target long enough so that the fire-control computer had completely 
charged the filters before commands were fed to the automatic pilot. 

This was accomplished by putting the following initial condition on the 
filter : 


e a (0) = E a (0) £ a (0) = o' 

e e (0) = E e (0) e e (°) = 


( 1 ) 


Because many combinations of autopilots, equations of motion, and 
aerodynamics were used, the different basic system configurations have 
been defined in table I. All configurations used will be referred to by 
the numbers given in table I. When nonlinear cross-product terms and 
aerodynamics that vary nonlinear ly with a are used, it will be stated, 
together with the configuration number to which the addition applies. 


RESULTS AND DISCUSSION 


All data-gathering runs were made with the target and guidance 
equipment included as part of the closed-loop system. Time histories of 
the airplane and fire-control variables were recorded from the time of 
radar lock-on to the firing time, t g = 0. 


Variables Used in the Analysis of Results 

Before discussing the results of this study, a brief discussion of 
some of the variables used in the analysis is in order. 

The direction cosines, which relate the interceptor body axes to 
space axes, determine the intercept or angular orientation. The direction 
cosines may be computed in two ways: by integrating the angular veloc- 

ities of the interceptor as was done in this study (eqs. (Al^)) or as 
trigonometric functions of the Euler angles (ref. 7)* When the latter 
method is used and the rotational order of the Euler angles is 0, 

$, the direction cosines l^ } m^, and n^, which relate the space Z-axis 

and the interceptor body axes, are given by 
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Zj = -sin 0 

rrij = cos 0 sin <£> > 

n^ = cos 0 cos $ 


(2) 


Thus, the direction cosine Z^ is always indicative of the pitch angle 
and for small values of 0, and n^ are indicative of the bank 

angle. However, n^ is considered more useful than m ^ as it is a 

cosine function and runs from -hi to -1 when the airplane rolls over onto 
its back. 


In addition to the direction cosines, the miss distances (eqs. (A 2 ) 
to (A5) ) or steering errors (eqs. (A6) and (A 7 )) are useful in analyzing 
the records. The miss distances are functions of the airplane motions 
and indicate, as a first approximation, the manner in which the inter- 
ceptor is maneuvering with respect to the target. Sometimes it is more 
convenient to use the steering errors, which are the miss distances 
divided by future range , to get a first approximation of the interceptor 
motion relative to the target. 

The usual airplane parameters u, v, w, p, q, and r and the 
control surface deflections 6 a , 5 e , and S r were also used to analyze 

the results. 


Basic Interceptor Flight Behavior With Linear Aerodynamics 

Flight maneuvers for various initial conditions .- Figures 6 and 7 
present the data obtained with basic configuration 1 for each initial 
condition. In these runs the airplane was represented by the complete 
six-degree -of-freedom equations of motion and linear aerodynamics. The 
low-gain longitudinal control was used. The direction cosines Z^ and 

n^, the normal acceleration n, and the steering errors, e a and e e , 

are shown in figure 6 and the linear and angular velocity time histories 
of the airplane are presented in figure 7 * As th e attack run always 
starts below the target, a climbing turn to the firing point is indicated. 
Figure 6 shows such a maneuver in which the interceptor had a normal 
acceleration of 5 S* In this case the normal acceleration and rolling 
commands were applied simultaneously . The major differences were the 
length of time the normal acceleration is applied, which depends on e e , 

-and the magnitude of the initial roll angle (see n^ time history), 

which depends primarily on the ratio G a | € e * As would be expected, the 

time of flight of the interceptor for each initial condition is different 
because of the differences in closing rate. 
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Characteristics of the motion .- In the normal acceleration time 
histories, the differences in the initial slopes are caused by differ- 
ences in the magnitude of the initial e e . The abrupt changes in slope 

that occur during the first 2 seconds of the normal acceleration records 
are caused by the g-limiter operation and rate-limiting effects in the 
longitudinal control system. As would be expected, these effects are 
clearly defined in the w and q time histories. In the case of q, 
for initial conditions II, IV, and V, the initial oscillations are partic- 
ularly violent because of the large initial e e value . 

When the interceptor has settled down to tracking, an oscillation 
occurs which is primarily evident in the p, v, and r records shown 
in figure 7 and the n^ record in figure 6. It was realized beforehand 

that some such oscillation was bound to occur because of the neglect of 
gravity in the roll command, since the neglect of gravity corresponds to 
commanding an uncoordinated maneuver. Consequently, when the interceptor 
banks it tends to drop off, that is, develop a lateral acceleration along 
the direction of the dipped wing. This lateral acceleration causes the 
interceptor to have a lateral velocity when e a approaches zero of such 

a sign as to cause the interceptor to "swoop" through the € a = 0 posi- 
tion. The interceptor then tends to reverse its bank to counteract the 
overshoot in e a , and some sort of oscillation may be expected to result 

in the final stages of the motion, involving primarily p, v, and r. 
Because of the gravity drop-off during the uncoordinated maneuver, the 
interceptor tends to come at the target from below, so that the final 
oscillation tends to be about the zero-bank condition. It may perhaps 
be roughly visualized as comparable to the swinging of a pendulum under 
the influence of gravity. 

A number of things may have a considerable effect on this oscillation, 
such as: various methods -of including gravity in the roll command, any 

parameters affecting the lateral accelerations while the airplane is 
rolling, the fact that the roll command tends to become indeterminate 
when both e a and e e approach zero, and the fact that the resolvers 

used in the simulator may tend to hunt when the roll command becomes 
indeterminate. In the present investigation there was not sufficient 
time available to evaluate the parameters affecting this oscillation. 
However, during the investigation of the effects of the nonlinear dynamic 
terms in the equations of motion of the airplane, it was found that the 
p term, which has an important effect on the lateral acceleration 
while rolling, also very strongly affects the final "swooping" oscil- 
lation. This phenomenon will be more fully discussed in a later section. 

A comparison of the data in figure 6 for each initial condition 
indicates that a qualitative similarity exists between initial condi- 
tions I and III and between initial conditions II, IV, and V. The time 
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histories of the airplane motion in figure 7 show that this similarity 
carries over to the airplane motions. Thus, two basic types of motion 
occurred, one where large changes in roll angles and small changes in 
flight-path angles were required to reduce the errors to zero (initial 
conditions I and III), and one where small changes in roll angles and 
large changes in flight-path angles were required to reduce the errors 
to zero (initial conditions II, IV, and V) . Because of this qualitative 
similarity, results presented in the following sections will be for two 
initial conditions, one for each basic type of motion. 

Comparison of yaw and sideslip dampers .- In the motions presented 
in figure 7, basic configuration 1 was used; thus, the yaw damper con- 
trolled the rudder. The basic configuration was changed to configu- 
ration 2 which replaced the yaw damper with the sideslip damper. Fig- 
ure 8 compares the v, r, and p motions of the airplane for these 
two rudder, controls. The (5 rudder control held the sideslip velocity 
to zero, except for a small excursion at the beginning of the run. The 
rolling velocity time history shows that p was not greatly affected 
by the change in rudder control. The yawing velocity is completely 
changed. The J3 rudder control yaws the airplane into the relative 
wind and in so doing maintains almost zero sideslip. The adverse yaw 
at the beginning of the run was eliminated and the peak yawing velocity 
almost doubled. The shape of yawing velocity time history when the side- 
slip damper was used resembled the shape of the sideslip time history 
obtained with the yaw damper. 

It is interesting to note that the change in rudder control had no 
effect on the terminal miss distances. 


Effects of Variable Aerodynamics 


Variable aerodynamics used .- As indicated in the section entitled 
"Components Simulated," the airplane stability derivatives were analoged 
to take into account variations with Mach number and angle of attack. 

The variation of the stability derivatives with Mach number was linear 
with the exception of which was quadratic. Equations (A20) and 

(A2l) describe the variation of the stability derivatives. The variations 


with a occur in C 




-'m* 


'ZpJ 


V and 


The angle -of -attack 


variation in 
coefficient. 


Cx is caused by the usual change in induced drag with lift 


In the case of C7 . a linear variation with angle of attack was 


included at all times in the basic C- 


This variation with angle of 


attack is a wing contribution. The nonlinear angle -of -attack variation 
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in Cj and the other stability derivatives are representative nonlinear 

variations of the stability parameters obtained from wind-tunnel tests 
of several configurations. The variations shown are not specifically 
for the configuration considered herein. Equations (A 2 l) show the form 
used to include these derivatives in the simulation. Since the nonlin- 
earities do not start at a = 0 , they were made functions of a - a cr 

where a cr is the angle of attack at which the nonlinearity starts. 
Provisions were made so that a cr could be varied. Figures 9(a) to 9(^) 
show the nonlinear variations employed for C m , C na , C^, and C 


for two different Mach numbers and a 


cr 


m> w npj> np 

= 3° and 6 °. Two forms of non- 


linearity in Cj were studied, a destabilizing one (fig. 9 (e)) and a 
stabilizing one (fig. 9 (f))* 


In order to define more clearly the effects of the nonlinear 
pitching moment, the time histories presented in this section are for 
basic configuration 4 , which uses the high-gain longitudinal control 
system. As the previously presented time histories are for the airplane 
with the low-gain longitudinal control system, it is necessary to present 
data for basic configuration 4 with constant aerodynamics to establish 
a standard so that the effect of the nonlinear pitching moment may be 
shown. 

Figure 10(a) (solid curves) presents the response of the airplane 
for basic configuration 4 with constant aerodynamics, initial condi- 
tion III. This record was cut off when a dropped below 3°, which 
corresponds to a time of 7 seconds . The addition of the nonlinear 
C m , a cr = 6°, to the simulation caused the airplane response to become 

unstable. (See dashed lines in fig. 10(a).) Because of this unstable 
condition, a cr is meaningless for this run and the time testing was 

cut off at 7 seconds . The saw-tooth character of the 6 e motion indi- 
cates that the instability is caused by a rate-limiting oscillation of 
the control surface . 


This same type of oscillation was encountered in the development 
of the roll control system and is discussed in reference 5 . Two methods 
by which the rate-limiting oscillation may be alleviated are to reduce 
the control surface deflection limit or add a q feedback to the con- 
trol system. Of these two, the q feedback is more effective and is 
preferred; however, it was somewhat more difficult to incorporate in the 
analog than a change of control-surface-deflection limit. Therefore, 
the latter was used in this study. 

Figure 10 (b) compares n, q, and 6 e for the 10° 6 e limit with 
and without the nonlinear C m , a cr = 6 °. An examination of the case 
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of linear C m shows that, when the 5 e limit is reduced, the ability 
of the airplane to develop normal acceleration is restricted and the aver- 
age normal acceleration is now about 4g with an overshoot to about 


on the initial transient. The control surface is now against the limit 
stops for about 8 seconds with a break at 1 second when the g-limiter 
kicked 5 e off the stop to limit normal acceleration. Because 6 e is 

now on the stop, the elevator is no longer available for damping. Thus, 
during the first 8 seconds an oscillation at the frequency of the airplane 
short-period longitudinal oscillation develops. The addition of the non- 
linear C m , which causes the airplane to become less stable as the angle 

of attack increases and corresponds to decreasing the static margin or 
increasing n|5 e , introduces changes in the longitudinal response. The 

nonlinear C m introduces a pitch-up, the effects of which are sensed 

.by the autopilot, and less control motion is required to obtain the 
desired normal acceleration. The elevator is now off the stop after the 
first three-fourths of a second and is able to damp the airplane short- 
period longitudinal oscillations. An examination of the normal accelera- 
tion presented in figure 10 (b) shows that, with the nonlinear C m , a 

higher average value of normal acceleration is attained. This causes 
the duration of the airplane normal acceleration to decrease. This higher 
normal acceleration was caused by the fact that the g-limiter operation 
depends on the n/& e response of the airplane for the high-gain control 

system and was not changed to account for the increasing n^ 5 0 introduced 
by the nonlinear C m . 


Results similar to the foregoing were obtained for all initial con- 
ditions with a cr = 3° or 6 °. Figure 11 compares the longitudinal air- 
plane response for the nonlinear pitching moment with a cr = 3 ° and 6 6 . 
When the nonlinearity was strong (a cr = 3°), the normal acceleration rose 

to higher values and the control surface displacement actually went to 
zero while the airplane was pulling maximum normal acceleration, because 
of the strong pitch-up characteristic. These same characteristics appear 
in the motions for a cr = 6 ° but the control-surface deflection never 

reaches zero when the nonlinear C m is effective. When the low-gain 
longitudinal control system was used, both values of <x cr gave results 

that paralleled those for the high-gain longitudinal control system. 
However, the effects were much less pronounced and initial conditions I 

and III were stable for = 20°. This occurred because the rate- 

e l 

limiting oscillation was more stable for the low-gain control system. 
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As pointed out in the foregoing discussion, the principal effect of 
the nonlinear pitching moment was to change the elevator motion in a man- 
ner that gives the desired command response . Because of the low time 
lags associated with the sensing instruments in the automatic control 
system, changes in the airplane motion, primarily the normal acceleration, 
are rapidly sensed and corrective signals generated. The high performance 
servomotors are able to move the control surface rapidly enough to com- 
pensate for these changes and maintain the desired response. 


Effect of nonlinear variations in the lateral stability derivatives . - 


The lateral stability derivatives, 


y v v “p 

tained nonlinearities that were functions of the angle of attack, failed 
to change the lateral response of the airplanes except when the addition 
included the C n nonlinearity with a™ = 3 ° 


C 1 , and 


Cn, which con- 
P 


cr 


The 


Cnp nonlinearity. 


like the other lateral nonlinearities, did not change the rolling veloc- 
ity response; however, the loss in directional stability caused the 
changes in the yawing and sideslip velocities shown in figure 12(a). 

The lack of change in the rolling velocity response is not surprising, 
as the guidance commands are calling for a specific rolling velocity and 
the changes introduced in the airplane response by the nonlinearities 
are picked up by the sensing instruments and a different control motion 
is ordered, taking into account these changes, that gives the commanded 
response. The changes introduced in the aileron response for nonlinear 


'V 


v 


and destabilizing are shown in figures 12(b) to 12(d). 


The other lateral nonlinearities produced no significant change in the 
aileron motion. The control motions for group additions of these various 
derivatives followed the trends indicated by the individual derivatives. 


Increasing a cr to 6° or a change to the low-gain flight-path con- 
trol system greatly minimized the effect of the nonlinearities in the 
lateral stability derivatives on the airplane and aileron motions . 


Effects of combined lateral and longitudinal nonlinear aerodynamics , 
Figure 13 presents the response of the airplane with all nonlinear aero- 
dynamics included in the simulation for <x cr = 3 ° and 6°. The attack 

run is defined by initial condition III and the simulation set up by 
basic configuration k. Changing a cr from 3 ° to 6° largely eliminates 

the effect of the nonlinear lateral stability derivatives. However, the 
nonlinear C m still produces a considerable effect. A comparison of 

the normal -acceleration curves presented in figures 11 and 13 for 
a cr = 3° shows that the combined lateral and longitudinal nonlinear 

stability derivatives, principally C and C , affect the airplane 


longitudinal motions to a greater extent than the nonlinear C m alone 
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The most significant differences are the negative normal acceleration 
developed in the combined case at approximately 6 seconds and the more 
oscillatory characteristics of the normal acceleration from 6 seconds 
to the end of the run. When a cr = 6° the nonlinearity is less extreme 


and the coupling of C m and 



nonlinearities is not as evident. 


A comparison of the sideslip velocities shown in figures 13(b) , all 
nonlinearities , and 12(a), lateral nonlinearities, for <x cr = 3° shows 


that the nonlinearity in C m affects the lateral response. The differ- 
ence between the solid and dashed curves of figure 12(a) is caused by 
the nonlinearity in C n . A comparison of the dashed curve of figure ll(a) 

P 

with the solid curve in figure 12(b) shows that, when the nonlinear C, 


m 


is present, large values of v are maintained longer. For example, at 
3.75 seconds, v = 190 feet per second for the nonlinear C n alone, 

H 

whereas it is 300 feet per second where the nonlinear C m is present in 
combination with the nonlinear C n p . Similar effects are produced in the 

yawing velocity. (Compare figs. ll(a) and 12(c) for a cr = 3°*) 


These changes in v and r are caused by the fact that the non- 
linear C m produces larger values of normal acceleration (about 8g com- 
pared with 5g when C m is linear) and the larger values of a make the 
Cnp nonlinearity more pronounced. 


Had the command g-limiter taken into account the changes in n^6 e 
of the airplane with the changing C m , no changes in v and r would 

have been noted as the normal acceleration and, consequently, the angle 
of attack would have been the same for the linear and nonlinear pitching 
moments. These data are discussed as ;they show some of the problems 
associated with a command g-limiter when the form of the nonlinearity 
in C m is not precisely known. 

Significance of nonlinear aerodynamic effects .- The primary result 
was that the guidance and control system was able to adjust itself to 
extreme variations in the aerodynamic characteristics of the interceptor 
and still cause it to perform the desired maneuvers and correct the lock- 
on guidance errors. The effects of loss of static longitudinal and 
directional stability ^represented by the C m and C n p nonlinear ities^ 

might have been more troublesome, however, had random external disturb- 
ances such as noise and gusts been considered. Unfortunately, lack of 
time available on the analog computer made it necessary to eliminate the 
investigation of random disturbance effects. 
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Moreover, even for the command response of the system there were 
significant secondary effects of the nonlinear aerodynamics on the control- 
system dynamics . For example, the nonlinearities had a considerable 
effect on the occurrence of rate-limiting oscillations in the control 
systems and on the operation of the command-type g-limiter assumed. 
Therefore, it is necessary to include nonlinear aerodynamic effects in 
simulating the command response of the automatic interceptor. 


The Dynamic Representation of the Interceptor 


Comparison of the complete and linear representations of the 
airplane . - The six-degree-of-freedom rigid body equations of motion, 
which will be referred to as complete equations, are: 


m(u + 

^0) 

+ fm [q.(Aw) - rvjj- = W 2 j + F x j 

m(v + 

ru Q - 

pw D ) + 

jm[r(£u) - p(zw)]| = Wmj + F Y 

m(w - 

^0) 

+ jmjpv 

- q.(z^u)]) = Wn 5 + F z 

J XP + 

[<* 

- i Y )qr 

= L 

i y 4 + 

[fr 

- x z)p r ] 

= M 

i z f + 

[( X Y 

- J x)p^ 

= N 


( 3 ) 


The linear equations used in this study were obtained by deleting the 
terms in brackets from the complete equations . 


As the complete equations accurately define rigid airplane motions 
and during its flight the rigid airplane actually solves the complete 
equations, differences between the interceptor motions for the complete 
and linear equations must be regarded as errors introduced by the use of 
a linear representation of the airplane. 

Figure Ik shows the motion of the center of gravity of the inter- 
ceptor and the bank angle for basic conf igurations 1 (complete equations) 
and 6 (linear equations). The first 3-5 seconds of an attack run defined 
by initial condition I are presented. In both cases the interceptor 
starts to bank toward the target, which is above and to the left of the 
t = 0 position of the interceptor. When complete equations are used to 
represent the airplane, the initial center-of-gravity motion is just what 
would be expected. It starts with a fall-off to the left due to gravity, 
but as normal acceleration develops the interceptor pulls upward to the 
left at the target. However, when linear equations are used to represent 
the airplane, the center of gravity moves to the right and swoops vio- 
lently upward as the airplane banks, before it starts to move left toward 
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the target. Thus, the use of a linear representation of the airplane 
introduces large errors in the center-of-gravity motion during the ini- 
tial hank and considerable errors in the bank angle. 

These differences in the motion indicate that large forces exist in 
the linear case that do not exist when the complete equations are used. 
These forces can be considered as spurious forces which are introduced 
because of the neglect of the nonlinear inertial terms enclosed in 
brackets in equations (3). Moreover, the initial upward swoop of the 
center of gravity caused by these apparent forces can be expected to 
greatly aggravate the previously discussed swooping oscillations, which 
occur in the later part of the attack run. An analysis of the cross- 
product terms which are neglected in the linear equations indicates that 
p(Aw), pv, q(ZV) are the important terms in the force equations, 
whereas pq and pr appear to be the important terms in the moment 
equations. Of these terms the lateral acceleration term p (At) and the 
normal acceleration term pv would have a direct effect on the initial 
center-of-gravity motion shown in figure 14 . The significance of these 
apparent or spurious forces can most easily be seen in the following 
manner . 

Consider the side-force equation, from equations (3), which is 
m(v + ru Q - pw Q ) 4- m£r(Au) - p(Aw)] = Wm^ + Fy 

and correctly states the side or Y-forces acting on a rigid body. When 
this equation is linearized under the assumption of small velocities and 
displacements the following equation results : 

m(v + ruQ - pw Q ) = Wm^ + Fy 

and correctly states the Y-forces acting on the rigid body when the 
assumptions of linearization are not violated. However, if the param- 
eters that are multiplied to give the neglected terms become large, and 
the magnitude of the neglected terms approaches or exceeds the magnitude 
of the retained terms, the linear equation is no longer a true statement 
of the Y-forces acting on the rigid body and the spurious or apparent 
forces occur. The nature of these spurious forces can be determined by 
changing the details of the method of linearization. Instead of deleting 
the cross-product terms, the linearization may also be accomplished by 
adding a force Fy* to the right-hand side of the general side -force 

equation that at all times is exactly equal to m(r Au - p(Aw)). Thus 
Fy * = m(r Au - p(Aw)) is the spurious force that occurs in the linear 

equation. 

In order to obtain some insight as to the magnitude of the terms 
being neglected when the cross products are deleted from the complete 
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equations pv , pq, pr, and q^Aw) were computed and are presented 

in figure 15 - As can be seen, very large forces are neglected at the 
beginning of the attack run. Because of the interaction of these terms, 
that is, p(Aw) affects w through pv which in turn affects v through 
p(Aw), the magnitudes of the cross-product terms under consideration for 
the complete equations are not the same as when linear equations are used. 
The cross-product terms which exactly represent the spurious forces in 
the linear case were computed and are presented in figure 15 . The most 
significant difference with regard to the incorrect initial center-of- 
gravity motion is that in the actual linearized run the p(Aw) term is 
considerably larger than the pv term. Therefore, the p(Aw) term 
appears to be of primary importance in accounting for the difficulties 
caused by the linear representation of the airplane. 

The interceptor system was next run on the simulator with the linear 
equations representing the airplane with the cross-product terms p(Aw), 
pv, pq, pr, and q(Aw) added in a systematic manner so as to determine 
their effect on the system and airplane responses . In addition, the 
effect of the other nonlinear terms was also checked during this part of 
the investigation. Figure 1 6 presents the corrections made by various 
cross-product terms to airplane motions, for initial condition I. A com- 
parison of figure 7 (a) with figure 1 6 shows that the terms p(Aw), pv, and 
pq definitely corrected errors in the airplane motion, and when these 
three cross-product terms are included the resulting airplane response 
is almost the same as it is for the complete equations. Inclusion of the 
p(Aw) term alone was sufficient to eliminate the gross errors, thus con- 
firming its relative importance as shown in figure 15. The term q(Aw) 
contributed a small but detectable loss in forward speed that was con- 
sidered unimportant. All other cross-product terms had little effect on 
the airplane motion. 

Figure 17 is included to show the effect of the nonlinear cross- 
product terms on the steering errors normal acceleration and direction 
cosine n^. There are considerable differences in these parameters 

during the attack run for the two representations of the airplane. At 
the firing time, t g = 0, a small azimuth error and large elevation errors 

exist for the linear representation, whereas both steering errors are 
almost zero when the complete equations represent the airplane . The 
steering errors that exist at t g = 0 for the linear airplane correspond 

to predicted terminal miss distances of J00 feet in azimuth and 500 feet 
in elevation. 

When the linear representation is used, the steering error in the 
latter part of the run remains relatively constant even though, as shown 
in figure 17(b), the airplane is pulling considerably more than lg during 
this part of the motion. Therefore, the aerodynamic load on the airplane 
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is apparently working against the spurious forces. In fact, a rough 
check of the average g-increment in the last 5 seconds showed that it 
was indeed very close to the average magnitude of the p(Aw) and pv 
values during this time . 

So far, the dynamic representation of the interceptor has been dis- 
cussed for the case of high rolling velocities, which corresponds to 
initial conditions I and III. In initial conditions II, IV, and V, the 
rolling velocity is much lower and since the important cross-coupling 
effects are functions of the rolling velocity one would expect different 
results for these initial conditions . An attack run for initial con- 
dition IV, which is a beam attack, will be used to illustrate the effect 
of reduced roll velocity on the dynamic representation of the interceptor. 

In this initial condition the maximum roll rate is 1.4 radians per second, 
compared with 3-5 radians per second for initial condition I. Figures 18(a) 
and 18(b) present the normal acceleration n, the direction cosine n^, 

and the steering errors e a and e e , and figures 18(c) and 18(d) compare 

the linear and angular motions for the complete and linear representations 
of the airplane when the attack run is defined by initial condition IV. 

Basic configuration 1 defines the simulator setup for the complete equa- 
tions and conf iguration 6 does the same for the linear equations . These 
time histories show that the roll cross -coupling effects are still pres- 
ent but the scale of the effect is greatly reduced. This is because 
the low rolling velocity causes the accelerations p(Aw), pv, and pq 
to have much smaller magnitudes and consequently the errors introduced 
by the omission of these terms is much smaller. In this attack run the 
increased loss in forward speed and increased Aw, compared with initial 
condition I (see figs. 7(a) and 7(d)), cause an increase in the impor- 
tance of two cross-product terms in the airplane motions. These terms 
are q(Aw) in the X-force equation and q(Au) in the Z-force equation. 

The q(Aw) term causes the velocity to fall off more rapidly and changes 
the shape of the Au curve (fig. 18(c)). The effect of q(Au) is to 
change the normal acceleration and w velocity over the midportion of 
the attack run. Figure 19 shows the effect of this term. These same 
effects were noted in initial conditions IV and V; however, in initial 
conditions I and III, the q(Au) term has negligible effect. 

The effects of the cross-product terms on the airplane . motion apply 
to all initial conditions . However, in initial condition V, the beam 
attack that degenerates into a tail chase, the final miss distances were 
small for both the linear and nonlinear representations of the airplane 
dynamics . It is felt that the combination of low rolling velocity, which 
reduces the coupling, and the long time of flight from lock-on to the 
firing point may be important factors in the difference in system effec- 
tiveness between initial condition V and the other Initial conditions. 
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Effect of dynamic representation on control systems .- An important 
interaction between the airplane and control system was found during 
this study. When the rate limit on control-surface deflection was reduced 
from 100 ° per second to 50 ° per second with basic configuration 6 (lin- 
ear equations) a rate-limiting oscillation occurred which caused the 
airplane to become unstable. For a similar reduction in the rate limit 
for basic conf iguration 1 (complete equations) the oscillation was greatly 
reduced and consequently the instability did not occur. The elevator 
and pitching velocities time histories for these two basic configurations 
are shown in figure 20. The important cross-product terms in this case 
were, as before, p(Aw), pv, and pq. However, the pv and pq terms 
are now of almost the same importance as the p(AwO term. 

The use of the 0 rudder control and/or the use of variable aero- 
dynamics did not change the results for the dynamic representation of 
the airplane. When nonlinear aerodynamics were included in the simu- 
lation, there was a definite tendency for rate-limiting oscillations to 
develop at 5 e ^ = 100° per second. Thus, as previously pointed out in 

the rate-limiting study, the pv and pq terms assume greater impor- 
tance than in the case of linear aerodynamics. 


Effects of Target Maneuvers 

A series of runs were made with a maneuvering target to determine 
the effect of target maneuver on the results pertaining to nonlinear 
aerodynamics and the dynamic representation of the airplane. In this 
study the target performed a ±2g vertical-plane maneuver. The target 
maneuver did not qualitatively change the results obtained for variable 
aerodynamics and the dynamic representation of the airplane with a 
straight-line target. 


CONCLUDING REMARKS 


The basic flight maneuver of the interceptor was a 5g climbing turn 
towards the target. In this maneuver, the pitch and roll commands were 
applied simultaneously . These attack-run results indicate that an auto- 
matic interceptor system that omits gravity corrections from the roll 
command can be made practical from a stability and control point of 
view, particularly if there is no overshoot in the longitudinal response, 
as was the case in the low-gain longitudinal control system. However, 
the large rolling motions in the final portion of the attack runs might 
present serious difficulties in connection with pilot comfort and radar 
tracking. More study is needed of the relationships between the command 
computer, automatic pilot, and guidance system. 
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When considered separately, the nonlinearities in the stability 
derivatives, except in the case of the nonlinear yawing moment due to 
sideslip C n produced no gross changes in the airplane motions. How- 
ever, the control-surface motion required to produce the commanded 
response was considerably altered by pitching-moment coefficient C m , 

C n , rolling moment due to sideslip , and damping-in-roll parameter 

P 3 


to a lesser extent . In the case of nonlinear 


C large changes 


occurred in the yawing and sideslip velocities, in addition to the changes 
introduced in the aileron motion. The combined nonlinearities in C m 

and C n showed definite coupling characteristics. Radar noise and 


atmospheric turbulence could have considerable effect on the airplane 
response when nonlinear aerodynamics are present. There is need for more 
research on this problem. 


The sideslip-angle rate rudder control proved effective in controlling 
the sideslip velocities encountered when linear aerodynamics was used in 
the simulation. The higher sideslip velocities introduced by the non- 
linear C n were also reduced by the sideslip damper. 


The study of the dynamic representation of the airplane showed that 
errors in airplane motion are introduced when a linear representation of 
the airplane is assumed. The addition of the cross-product terms p(Aw), 
pv, pq, and pr to the linear equations corrected these errors. Of 
these, the p (Aw) term was most important. When the magnitude and time 
of application of normal acceleration increased, the cross-product terms 
q(Au) and q(Aw) became more important, indicating that for a maneuvering 
target these terms may become important. 

A vertical -plane target maneuver of ±2g did not change the results 
for variable aerodynamics and the dynamic representation of the airplane 
obtained with a straight-line target. 


The results presented in this report are for one type of automatic 
interceptor system and are limited in Mach number and initial condition. 
Thus, these results are of limited generality, and extensions to airplanes 
in other interceptor systems should be made with caution. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., September 25, 195&- 
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APPENDIX A 


DISCUSSION OF THE ATTACK PHASE OF THE AUTOMATIC INTERCEPTOR 
SYSTEM SIMULATED ON THE TYPHOON COMPUTER 


Reference 8 is a comprehensive review of the airborne fire-control 
problem. The general theories of airborne fire control are developed 
in this reference. This appendix is limited to a discussion of the 
specific interceptor system set up for the Typhoon computer. 

Figure 21 is a detailed block diagram for the simulation of this 
interceptor system. Information and information flow routes and the 
relationship of the equations to each other and as part of the inter- 
ceptor system are indicated by this diagram. The numbers within the 
blocks refer to equations in this appendix. 

The target (figs. 1 and 2) was programmed to fly a straight-line 
course at an altitude of 50,000 feet or perform a ±2g vertical -plane 
maneuver that starts at the same altitude. When the maneuver was used, 
it was started at radar lock-on. The target flight was along the earth 
Y-axis in a positive direction or confined to the YZ-earth plane during 
the maneuver. The following equations were used to simulate the target 


where 



l 


1 


l 


3 


T - 1 


m^ m 2 


n l 


n 2 



V T = 


^ 

i (0) + 


j Vqi cos 7ijt 


— ^ 

k Viji sin ?ijit 
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and 
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7 t = “vj- 
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where 

= 0 for a nonmaneuvering target 

rvp = +1 or +2 for an upward-maneuvering target 

nrp = -1 or -2 for a downward-maneuvering target 

The coordinate transformation T - 1 as indicated consists of the 
direction cosine matrix, the direction cosines being obtained from equa- 
tion (Al4) . Equation (Al) converts the target velocities from earth 
axes to the interceptor principal body axes. In this problem Vip was 

a constant. 


The radar was assumed to be dynamically perfect and therefore could 
be simulated by using the geometry and kinematic relations between the 
interceptor and target. The range rate is given by 



u 

v S 
w 


(A2) 


where u, v, and w are the components of the interceptor velocity 
vector in principal body axes. Then the components of R are given by 



and the radar gimbal angles 


Q e , and R are given by 
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R x sin 0 a - Ry cos 9 a = 0 


R x cos 0 Q + R v sin 9 


a ^ 


- ^x 2 + 


/ 


H x 2 + Ry 2 sin 0 e + R z cos 0 0 = 0 


(A4) 


v / % i T 


Ry 2 cos 0 e - R z sin 0 e = R 


The fire-control computer was assumed a director -type computer, which 
uses the range and rate information to solve the equations of a lead 
collision rocket-firing course. In setting up the solution of these 
equations, it was assumed that no miss existed along the radar line sight, 
which is equivalent to assuming a two-component miss-distance vector in 
radar axes. The equations solved by the fire-control computer are 


R 


t g + t 


V T 


R x - . 

X t~ + T 


cos 0 a cos 0 e + 


Rv - 


PV 
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Rtj 


sin 0 a cos 0 e 


Ry. - 




t g + T . 


sin 0 e = 0 


t g + T t g + T 


R- 


PV 




M a 


Rr 


t g + T 


t g + T 


R Z - 


aV m T 


t g + T 


(A5) 


These equations are written in a mixed axes system. The first equation 
is written in radar axes and the last two in interceptor axes. This was 

M M 0 

done so that the output of the guidance computer, - — and 


t g + t 


t g + t 


could be used to compute automatic pilot commands without making a trans- 
formation from radar axes to interceptor axes. It should also be noted 
that these equations were divided by tg + T. This was done to improve 

scale factors and voltages in the computer simulation. 
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The 

steering 


fire -control- computer output was converted to unfiltered 
errors by the following equations: 


E Mtg + T ) 

v""t 
V f + - ■ w- 

1 t„ + T 


E e - 
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V T 

* w 
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'm 


f t g + T 


(a6) 


These unfiltered steering errors were then passed through the filter and 
corrected for cross roll. These modified steering errors are given by 


e a " 1 A fl b ( E a - T « pe e) 
e * “ + T ^ 


(A7) 


The that comes from equation (A 7 ) is the unlimited vertical flight- 

path change command. The e e is passed through a g- limiter that 

restricts the command to a predetermined level. The equation of the 
g-limiter is 


< 6 P 

e min e e max 


(A8) 


and the output is identified as e e which is the limited value of € e . 

The minus sign is attached to € P . because it was computed as a posi- 

min 


tive number. 


The limits used in equation (a 8) are independent of the airplane 
and automatic pilot characteristics when the low-gain flight-path control 
system is used, whereas the limit is dependent on the airplane and auto- 
pilot parameters when the high-gain flight-path control system is used. 
These normal acceleration limits were computed by the following formulas 
for the flight-path control system: 
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For the low-gain system 

e max 

For the high- gain system 

_ HlS^ + k 12 + 3-96) _ n£g^C6 + K 12 + 3-96) , . 

e max VpK^Kg e min VpK^Kg v 

The output of the g-limiter e e (see eq. (A 8 )) is the command, to the 
flight-path control system. 

The roll command was computed as 


n L g 

V* 


-mm 


n Lg 

V F^ 


(A9a) 




tan 


-1 


e a 


(A10) 


As the outer feedback loop for this control system is through the 
radar and computer, the inputs to the automatic pilot, e e and $-p, are 

used to command a flight-path angular velocity and rolling velocity. 


Gravity was neglected in computing the roll and flight-path commands 
in this study. This was done because it was one of the objectives of 
this study to determine if a successful system could be set up neglecting 
these gravity corrections. 


The automatic pilot used in this study is represented by the 
following equations: 

■ ( 1 + T SE D )5 e - A + ^)( K 4«e - r) + K 12l + k 20 


(! + T SA D ) 5 a " -^c + + 


(All) 


The rudder commands, which were regulatory in nature, were given by 
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or 


(1 + t sr D) 5 r = K u r ' 
( 1 + T SR D ) 5 r = " K llP v 


(A12) 


The block diagrams of these control systems are presented in figures 3 
and 4. 


As indicated in the equations and system constants given at the end 
of this appendix, the roll control system is characterized by a high for- 
ward loop gain, a first-order servo, and a rolling velocity and rolling 
acceleration feedback. The rolling velocity feedback is for damping and 
to close the inner loop of the control system. The rolling acceleration 
feedback is included to relieve rate-limiting oscillations (ref. 3 ) • 

The flight-path control system has a first-order servo and a 7 , flight- 
path rate, feedback to close the inner control loop and a 0 , pitching 
velocity, feedback is included to supply damping. A low-gain system 

y 

that gives — = Q.4 was used for nonmaneuvering target and a high-gain 
e e 

y 

system with — =1.4 for the maneuvering target as well as for the 
G e 

case of the nonmaneuvering target. The rudder is controlled by either 
an r, yaw-rate, or (3, sideslip-angle rate, feedback. The yaw-rate feed- 
back is used to damp the Dutch roll mode and is not effective in con- 
trolling the sideslip. Therefore, the airplane turn is poorly coordi- 
nated when the r feedback is used. The 0 feedback provides Dutch 
roll damping and, in addition, improves the turn coordination. 

The servomotors were simulated as first order with a time lag of 
0.03 second. Provisions were incorporated in the servomotor analog so 
that the maximum rate and displacement of the control surface could be 
limited . 

With the high gains used in these control loops, there is a possi- 
bility that second-order servomotors might produce an unstable control 
system for the same conditions that a control system with first-order 
servomotors would indicate stable conditions. Accordingly, an analysis 
was made using a second-order servomotor which would have an equivalent 
time lag of 0.03 second to first, order. The damping ratio assumed 
was 0.35, and the natural frequency was 3*72 cycles per second. The 
results of this analysis indicated that the rolling velocity and flight- 
path rate -control systems as used in this interceptor study would be 
stable with second -order servos . 
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The use of a flight-path control system is inconsistent with the 
omission of gravity corrections as the most practical way to measure 7 

is 7 = “ gHj)- As originally planned, the longitudinal control 

system was a 0, pitching -velocity, control system. The 7 control 
system was substituted for the 9 control system when the vertical-plane 
studies of the attack problem indicated that a 7 control system would 
be more satisfactory and easier to develop than the 0 control system. 
Unpublished results from vertical -plane studies, made during and since 
the Typhoon work, have indicated that a 0 control system as satisfactory 
as the 7 control system should be obtainable. 

Since 7 could be obtained on the simulator without a vertical 
reference and because it was one of the objectives of this study to find 
out if a system without a vertical reference was possible, such a device 
was not included as part of the simulation. 

For the Typhoon study the interceptor was represented by the six- 
degree -of -freedom rigid body equations of motion. Engine momentum terms 
were neglected as ramjet power was assumed for the attack run. The equa- 
tions of motion were referenced to principal body axes, thereby making 
the product of inertias zero. Thus, the equations of motion become 


m(u 

+ l w o) 

+ jm fq (/w) - rv]j 

m(v 

+ ru D ■ 

■ Pw c ) + |m[r(Zhi) 

m(w 

- <JP 0 ) H 

- jm[pv - q(Au)]J 

! X P 

+ [Uz 

1! 

V 

1 

I Y q 

+ [(% 

•- *z)p^] = M 


♦[(* 

- I x )pq] = N 

where the terms 

in the 

brackets are the 


(A13) 


that are deleted to obtain the linear equations . These equations were 
analoged so that all or any one or combination of the nonlinear could 
be deleted from the simulation. The linear equations used in this study 
differ slightly from the classical linear equations in that the direction 
cosine computation was not linearized and, therefore, the gravity terms 
in the force equations contain no linearizing assumptions as they do in 
the classical linear stability equations . 
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The direction cosines were computed by the following equations: 


l± = l± + (mqr - nqq)dt 

0 J 0 

m i = m i Q + J ( n i p " z i r ) dt [ (i = 2 > 3) (Al4) 

n i = n i 0 + f Q ( l l<l ~ m iP) dt 

The aerodynamic forces and moments are Fy, Fy, and Fy and L, M, and 

N of equations (A13) and consist of the control-surface forces and 
moments, the damping forces and moments, and the static forces and moments. 

The control-surface forces and moments were computed by 

q s M 2 S 


( 4 e - c v* 


m 


(5) 


m 


q Q M 2 S 


- C Z 5 /e — 


(A15) 


fe) 

W 


= f C 1 . & a + C Z* 5 r 


V* 


2 sb 


&a> ^r 


5 r 


N 


q^M 2 S5 

- V“ 

q s M 2 Sc 


= C„ 5 


'8r 


n c '■'r t 
&r 2 


(Al6) 


The r ema ini n g forces and moments which are caused by the airplane motions 
are computed by 
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X _ X Q [ C x q s M2s 


ni m 


m 


M^S 


Y 

m m 


z c % m2s 

m " m 


(A17) 


H 


2V S ' 2V S 


bCi 

r r + MC 7 3 

t-D 




M 


N _ | bCn P 


mq 


1 


q. + 


MCjji 


j 


MSc 

qs IT 


b 

2 


I Z I 2V s 


p + ^ + M S p ] qs i 


(A18) 


and for equations (A13) 


^ TT ^ 

fx 

m 


m 


= (A15) + (A17) 


and 


m 

^ s 




ix 

M_ 

lY 

N_ 

l>J 


= (ai6 ) + (Ai8) 


The problem required that the angle of attack a, the sideslip 
angle $, and Mach number be supplied, as well as values for the stability 
derivatives . The following equations were used to compute <x, 3, and 

Mach nu mb er, as well as V, (3, and 7 which were used in the guidance 
and automatic pilot: 
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7 = 


F z + e n 3 

mV 

J 


(A 19 ) 


The aerodynamic stability derivatives, except c x> and the control 

surface effectiveness were determined by analytical investigation and 
the results were published in reference 9 * These stability derivatives 
are the linear derivatives and are functions of Mach number . The control 
surface effectiveness and the drag coefficient C-^- were obtained from 


unpublished wind-tunnel tests. Thus, the drag coefficient is a nonlinear 
function of angle of attack, as well as a function of Mach number . In 
addition, the basic CL is composed of two parts, a contribution from 


the wing and a contribution of the fuselage and tail. The contribution 
to C^p from the wing is a function of angle of attack. Therefore, the 


basic C 7 


is a function of angle of attack, as well as of Mach number. 


The pitching-moment coefficient C m and the stability derivatives 

C ? , CL , CL , and CL were programmed with a nonlinear tail that is 

P P 

a function of the parameter a - a cr - Here a is the angle of attack 
from zero lift and a cr is the angle of attack at which the nonlinear 
tail starts. Provisions were so that a cr could be varied; two values 
were used in this study, 3° and 6°. 

All stability derivatives used were plotted as functions of the Mach 
number and angle of attack. These curves were then fitted by least 
squares, to obtain a polynomial function of Mach number and angle of 
attack that did not vary more than ±10 percent from the original curve . 
This procedure led to the following equations, which were used on the 
Typhoon computer to make a continuous determination of the stability 
derivatives as the Mach number and angle of attack changed 



3 6 


NACA EM L56JO8 


M*C X = (0.4437 - 0.0335M)a' 
°e 

M*C y = 0.272 
5 r 

M*C Z = -0.1514 

► 

M*C 2 = 0.0537 
6 r 

= -0.8438 
M*C n& = -0.2297 


and. 

M*C X = -0.2425 - 0.1325M* + (47.105 - 38.364m* + 7.38M* 2 )a 2 


M*C Y = 0.7644 - 0.284M* 
3 


M*C Z = -5.15a 


M*C„ = -0.4437 + 0.0535M* 

C Zp = -0.404 + 0.0975M* ± 50.7(a - a c )2>i(a - a cr ) 


c Zr = 0.2689 - 0.0565M* 


M*C Z(j = -0.1957 + 8.10(a - a c )%i(a - a cr ) + 


0 4 = -6.320 + 3 .74M* (M* < I.7988) 


0 4 = 0.913 - 0.281M* (M* > I.7988) 


(A20) 


(A21a) 

(A21b) 

(A21c) 

(A21d) 

(A21e) 

(A21f) 

(A21g) 

(A21h) 

(A21i) 


= -3-959 + 0.535M* 


(A21 j) 
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M*C m = -1.369a + 1^31 .91 (a - a cr )2 - 152. 4 (a - a cr )3 

- a . cr J 

(A21k) 

c np = - 0.1864 + O. 1625 M - 0 . 0281 M * 2 - 1 . 197 (cl - a cr ' 

) 2 ^( a - a cr) 

(A 21 Z) 

c n r = -0-955 + 0.18M* 


(A21m) 

M*C np = 0.528 - 13. 97 (a - a cr )2^ a - 

a cr) 

(A21n) 


where 


p ^cx — cCq-joJ " 0 ,a ^ ^cr 

p ^00 “ a cr ^ = > co cr 

These equations for the aerodynamic forces and moments contain all of the 
necessary information to obtain continuous aerodynamic information as the 
interceptor angle of attack and Mach number change . These equations were 
set up on servo multipliers and arranged so that any of the following 
combinations could be used: 

(a) constant aerodynamics 

(b) aerodynamics varying with Mach number 

(c) aerodynamics constant with Mach number and varying with a - a cr 

(d) aerodynamics varying with Mach number and a - a cr 


The parts of these equations that vary with a - a cr were zero for 
a < a cr and for a > a cr the computed value was used. The equations 
were so analoged that the value of a cr could be changed manually. In 
addition, the a ?■ term in and the a term in C ^ could be removed 

if desired. The stability derivatives C^, C n ^, C^, and were 

neglected in this study. 


In order to include the nonlinear aerodynamics, variations of air 
density with altitude had to be neglected. Thus, even though the 
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interceptor changes its altitude "by several thousands of fe.et, the air 
density remains constant at the value for 50,000 feet, the initial alti- 
tude of the target. 


The physical constants used in the problem and airplane character- 
istics are given in the following tables: 


Constants : 

g, ft/sec 2 

¥, lb 

m, slugs 

I]_, slugs -ft 2 

I2, slugs-ft 2 

Ij, slugs-ft 2 

S, sq. ft 

b, ft 

C, ft 

Vgi, ft/sec 

Vp, ft/sec 

V s , ft/sec 

V I Q > ^/sec 

qg (at H = 50,000), lb/ft 2 

r, sec ..... 

t j sec •>••••••• 

t sa> t se> t sr > sec • • • • 

<x cr , deg 

^and^ 

m i z 


32.17 

25, 000 
777 
15,040 

274,800 

287,000 

401 

35-81 

15 

1359.5 

2,000 

971.1 
2,140 

170.2 

1-5 

0.6 

0.03 

3 or 6 


Values required for 

£ = 0 and -SL = 0 at t = 0 
m 


Flight -path autopilot: 


*6 

k 8 

K 12 

%) 


Low gain 
0.4 

14 

20 

2 


High gain 
2 

14 

0 

1 


Adjustable constant, 0 to 0.2 
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Roll autopilot: 

%••••'• 

k 9 

K 10 • • • • 


... 8 
... 5 

0.1 to 0.2 


Rudder autopilot: 

K n 

K n 


Configuration 1 

1.8 

0 


Configuration 2 
0 


- 1.8 


Limiters : 

5 a , 5 e , 5 r , deg limited at ±10 or ±20; normal value ±20 

S a , 6 e , 5 r , °/sec limited at ±20, ±50, 

±100, ±143; normal value, ±100 


g-limiter : 


e e 

c max 


Low-gain autopilot 
(*) 
4g_ 


Vi o K4 


High-gain autopilot 

(*) 

4g K 6 + Kj_g + 5 .96 

Vi 0 K4 K6 


£e min 


2g 4g 1% + K lg + 3.96 

V I Q K 4 V I q k 4 *6 



V 


I* 


is initial value of 


ho 
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APPENDIX B 

THE EFFECT OF AUTOMATIC -PILOT FEEDBACKS 
ON AIRPLANE CHARACTERISTICS 


The automatic pilots used in the Typhoon study contain both angular 
rate and acceleration feedbacks. The use of these feedbacks effectively 
changes the inertia and aerodynamic damping characteristics of the air- 
plane* The following table shows the effect of the rolling velocity and 
rolling acceleration, yawing velocity, and pitching velocity on the airplane 



Basic 

airplane 

Change produced 
by autopilot 

Feedback 

producing 

change 

lx 

15,040 

81,400 

P 

p ! 

-0.1895 

-24.6 

P 

C l 

°r 

0.1446 

5.24 

r 

% 

-0.559 

-22.4 

r 

C™ 

“q. 

-2.782 

-218 (low-gain longi- 
tudinal autopilot) 
-109 (high-gain longi- 
tudinal autopilot) 

q. 

q. 
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TABLE I 

BASIC SYSTEM CONFIGURATIONS FOR DATA GATHERING RUNS 


^Target and guidance included in all basic configurations; control 
surface rate and displacement limits are 20° per second 
and 100° per second, respectively, unless otherwise notedj 


Configuration 

Airplane 

equations 

Lateral control 

Longitudinal 

control 

Aerodynamics 

1 

Complete 

Roll pilot with 
yaw damper 

Low gain 

Constant 

2 

Complete 

Roll pilot with 
sideslip damper 

Low gain 

Constant 

3 

Complete 

Roll pilot with 
yaw damper 

Low gain 

Varies with 
Mach number 

4 

Complete 

Roll pilot with 
yaw damper 

High gain 

Varies with 
Mach number 

5 

Complete 

Roll pilot with 
sideslip damper 

High gain 

Varies with 
Mach number 

6 

Linear 

Roll pilot with 
yaw damper 

Low gain 

Constant 
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Figure !•- Information flow diagram for the attack phase of the automatic 

interceptor. 
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Figure 3.- Block diagram of lateral control system. 
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(a) Sketch showing the parameters used in defining 
and a table of parameter values for the initial 


initial conditions 
conditions studies . 


Figure 5.- Initial conditions for the investigation of the attack phase 
of the automatically controlled interceptor. 
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t , sec 

(a) Initial condition I. 

Figure 6.- Direction cosines, normal acceleration, and steering errors 
for the attack i*un. Interceptor represented by complete equations 
with constant aerodynamics. Basic configuration 1 . 






p,!0q,)0r, rad/sec 
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Figure 8.- Comparison of lateral interceptor response for yaw and side- 
slip dampers. Interceptor represented by complete equations and con- 
stant aerodynamics. Initial condition I; basic configurations 1 and 2 
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(b) 8gj — 10° j — 


Figure 10.- Concluded. 
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(d) Effect of nonlinear , destabilizing, on the aileron response 

ir 


Figure 12.- Concluded. 
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(b) Linear velocities. 
Figure 13 . - Continued. 
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(c) Angular velocities. 
Figure 13 •- Continued. 
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(b) Normal acceleration and direction cosine n^. 


Figure 17 .- Concluded. 
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(a) Normal acceleration and steering errors. 

Figure 18. - Comparison of the interceptor response for the complete and 
linear equations of motion. Low rolling velocity case; initial con- 
dition TV; constant aerodynamics; basic configurations 1 and 6. 
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(d) Angular velocities. 
Figure 18.- Concluded. 
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